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We study the measurement for the unambiguous discrimination of two mixed quantum states 
that are described by density operators pi and p2 of rank d, the supports of which jointly span a 
2d-dimensional Hilbert space. Based on two conditions for the optimum measurement operators, 
and on a canonical representation for the density operators of the states, two equations are derived 
that allow the explicit construction of the optimum measurement, provided that the expression for 
the fidelity of the states has a specific simple form. For this case the problem is mathematically 
equivalent to distinguishing pairs of pure states, even when the density operators are not diagonal in 
the canonical representation. The equations are applied to the optimum unambiguous discrimination 
of two mixed states that are similar states, given by p2 = UpiW , and that belong to the class where 
the unitary operator U can be decomposed into multiple rotations in the d mutually orthogonal 
two-dimensional subspaces determined by the canonical representation. 
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I. INTRODUCTION 

The discrimination of nonorthogonal quantum states 
is of fundamental interest for many problems connected 
with quantum communication and quantum information 
and has consequently attracted a great deal of attention. 
An overwiew of the theoretical aspects of quantum state 
discrimination is given in recent review articles [3, Q]. 
The standard problem is the following: We assume that 
a quantum system is prepared in a certain state that is 
drawn with known prior probability from a finite set of 
known possible states, and we want to find the best mea- 
surement for determining the actual state of the system. 
When the given states are nonorthogonal, they cannot be 
discriminated perfectly, and therefore various measure- 
ment strategies have been developed that are optimized 
with respect to different criteria. The most prominent of 
these schemes are discrimination with minimum error Q 
on the one hand, and optimum unambiguous discrimina- 
tion 0, [fl on the other hand, and very recently also the 
strategy of discrimination with maximum confidence has 
been introduced 

In a measurement for unambiguous discrimination er- 
rors are not allowed, at the expense of admitting incon- 
clusive results, where the measurement fails to give a 
definite answer. In this paper we restrict ourselves to 
considering only two given states that in the most gen- 
eral case are mixed states. Clearly, error-free discrim- 
ination of mixed states is only possible when the sup- 
ports of the states are not identical. Note that the sup- 
port of a quantum state is the Hilbert space spanned by 
those eigenvectors of its density operator that belong to 
nonzero eigenvalues, and the rank of a state is the dimen- 
sion of its support. The optimum error- free measurement 
we are trying to find is the measurement that minimizes 
the average probability of getting an inconclusive result, 
or in other words, the average failure probability, where 
the prior probabilities for the occurrence of the differ- 



ent possible states are taken into account. We mention 
that recently unambiguous discrimination was also in- 
vestigated without considering these prior probabilities, 
by requiring that in the best measurement the largest 
state-selective failure probability for any of the incoming 
states be as small as possible [7|- Here we stick to the 
traditional way of defining optimality for unambiguous 
discrimination by requiring that the average overall fail- 
ure probability of the discriminating measurement be as 
small possible. 

While the optimum measurement for the unambigu- 
ous discrimination of two pure states was found already 
a long time ago 0, unambiguous discrimination in- 
volving mixed states, or sets of pure states, equivalently, 
became an object of research only more recently [8-20]. 
So far a general analytical solution for the optimum mea- 
surement that unambiguously distinguishes between two 
arbitrary mixed states does not exist yet, but a number 
of general results have been obtained. Several necessary 
and sufficient conditions for the optimum measurement 
have been derived [1, and it has been shown that 
the solution can be found in an efficient way using the 
method of semi-definite programming [9|1. Moreover, re- 
duction theorems have been developed [lOl [20| that can 
simplify the discrimination problems. 

Lower bounds for the failure probability [IlL E3, Eil 
as well as the conditions for saturating the bounds 
HU EH Gil have been also studied. It has been estab- 
lished that the minimum failure probability in the unam- 
biguous strategy is at least twice as large as the minimum 
probability to get a wrong result when errors are allowed 
to occur [151 ] . As a consequence of the reduction theorems 
[8] it follows that the overall lower bound of the failure 
probability, proportional to the fidelity [Tl| , can only be 
saturated when the dimension of the joint Hilbert space 
spanned by the supports of the two states is equal to 
the sum of their ranks. Even for this case the saturation 
of the bound depends sensitively on the structure of the 
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density operators and on their prior probabilities, and it 
is expected [H, EH that in general the fidelity bound can 
be reached only in a very limited range of all parameters. 

In a few special cases a complete analytical solution 
for the optimum measurement, valid for arbitrary prior 
probabilities of the two given states, has been derived. 
These cases are the unambiguous discrimination of 

i) a pure state and an arbitrary mixed state, known as 
quantum state filtering [l6|, [TtJ , 

ii) two density operators of rank d in a (d+l)-dimensional 
joint Hilbert space [ll|, and 

iii) two density operators of rank d that are simultane- 
ously diagonal in the canonical basis [ll| that separates 
the 2d-dimensional joint Hilbert space into d mutually 
orthogonal two-dimensional subspaces [l8[ . 

The solution derived in Ref. [l8| includes the unambigu- 
ous discrimination of two uniformly mixed states, which 
is equivalent to the discrimination of the two subspaces 
spanned by their supports. Moreover, the case iii) also 
applies to the comparison of two given pur e states (2~0 | 
having arbitrary prior probabilities l22j| and to the 
programmable discrimination [23| that distinguishes be- 
tween two pure states when one [24[ or both of them 
[23l H3, are unknown and the discrimination is per- 
formed with the help of reference copies. 

In the three cases listed above the optimum measure- 
ment can be constructed from the solutions for discrim- 
inating pairs of pure states. When such a reduction to 
pure-state discrimination problems is not possible, two 
classes of analytical solutions are known. First, general 
expressions for the optimum measurement operators have 
been derived which hold in the special case that for the 
given prior probabilities of the states the lower bound of 
the failure probability is saturated [3] . Second, the opti- 
mum measurement has been determined for two equally 
probable geometrically uniform states of rank two, de- 
scribed by density operators p% and p 2 with p 2 = Up\U, 
where U 2 — I with / being the identity [2fJ. 

In the present paper we extend the above list of cases 
where a complete solution, valid for arbitrary prior prob- 
abilities of the states, can be obtained. We show that 
when a certain condition for the fidelity of the two states 
is fulfilled, the solution of the case iii) can be used to 
solve the optimization problem also for states that are 
not diagonal when represented with the help of the re- 
spective sets of orthonormal canonical basis states. In 
particular, the required condition for the fidelity is found 
to be satisfied for two mixed states that are similar, given 
by P2 — UpxW where WU = I, and that in addition be- 
long to the class where the unitary operator U can be de- 
composed into multiple rotations in the two-dimensional 
subspaces determined by the canonical representation of 
the mixed states. Interestingly, by extending the quan- 
tum key distribution protocol based on two nonorthog- 
onal pure states [2(| to the case of two mixed states, it 
has been found a decade ago [27[ that secure commu- 
nication is only possible in this protocol when the two 
mixed states are connected by a rotation operator with 



a nonorthogonal angle and belong to the special class of 
states considered in this paper. 

The paper is organized as follows: In Sec. II we review 
earlier results that are needed for the present investiga- 
tion, and we derive our basic equations. The optimum 
measurement for the unambiguous discrimination of the 
special class of states considered in this paper is obtained 
in Sec. Ill, and the paper is concluded in Sec. IV. 

II. GENERAL THEORY 

A. Conditions for the lower bound of the failure 
probability 

We start with a brief summary of the basic theoretical 
concepts and results that are needed for our further treat- 
ment. Any measurement for distinguishing two quantum 
states, characterized by the density operators p\ and p 2 , 
can be formally described by three positive detection op- 
erators obeying the equation 

Il + Il 1 +Il 2 = I, (2.1) 

where / is the identity. These detection operators are 
defined in such a way that Tr(/oII k ) with k = 1, 2 is the 
probability that a system prepared in a state p is inferred 
to be in the state pk, while Tr (pHo) is the probability 
that the measurement fails to give a definite answer. The 
measurement is a von Neumann measurement when all 
detection operators are composed of projectors, other- 
wise it is a generalized measurement based on a positive 
operator- valued measure (POVM). From the detection 
operators life schemes for realizing the measurement can 
be obtained using standard methods [28] . For the results 
of the measurement to be unambiguous, errors are not al- 
lowed to occur so that there is never a misidentification 
of any of the states. This leads to the requirement 

Pl n 2 = p 2 U 1 = (2.2) 

[H Q, which means that Tr(p k II ) = 1 - Tr(p k II k ) for 
k = 1,2. When we denote the prior probabilities for the 
occurrence of the two states by 771 and rj 2 , respectively, 
with 771 + 772 = 1, the total failure probability of the 
measurement, Q, is given by 

Q = 7yiTr(piilo) + r/2Tr(p 2 no) 

= 1 - TnTrO^rii) - r72Tr(p 2 n 2 ). (2.3) 

From the relation between the arithmetic and the ge- 
ometric mean we get Q > 2-\/?7i?7 2 Tr(/9ino)Tr(/9 2 no), 
and because of the C auchy- Schwar z- inequality this yields 
Q > 2 v fffT^Ma,x v \Ti(Vy/pinoy/pZ)\ 0, where V de- 
scribes an arbitrary unitary transformation. The failure 
probability takes its absolute minimum when the equal- 
ity signs hold in these two relations, which is true if and 
only if both the equations 

^TrO^Ho) = 77 2 Tr(p 2 II ) (2.4) 
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and V^piy/Tlo ~ ^/ply/^o are fulfilled. After multiply- 
ing the second relation with its Hermitean conjugate, the 
two conditions for equality can be combined to yield [l3[ 

771 x/n^pi Vno = ^Vn^Vrio- (2.5) 

Substituting IIo = 7 — 111—112 into the inequality for the 
failure probability Q, given above, we arrive at 



where 



Q > Zy/mm F(pi,p 2 ), 



F = Tr [(V^ PiVj^) 1/2 } = Tr| VpTv 7 ^ 



(2.6) 



(2.7) 



is the fidelity [23]. From Eqs. (J23|) and JZ3]) we conclude 
that the lower bound of the failure probability, propor- 
tional to the fidelity of the states, is obtained if and only 
if 77iTr(piIlo) = 772Tr(/9 2 n n ) = y/rjTrj^F. This is equiva- 
lent to the two conditions [13| 

Tr(p 1 H 1 )-l + x [^F(p 1 ,p 2 ) = 0, (2.8) 

V m 

Tr( /02 n 2 )-l + x f^F(px,p 2 ) = (2.9) 

V V2 

that are the basic equations for our further treatment. 
Whenever we can find detection operators 111 and II2 
satisfying Eqs. (gUJ) and ((2U) while IIq = I - Hi - II 2 
is also a detection operator, i. e. a positive operator 
with eigenvalues between and 1, then we are sure that 
these operators determine the optimum measurement for 
unambiguously discriminating the states, since they yield 
the lower bound of the failure probability, proportional 
to the fidelity. In the optimum measurement the lower 
bound can only be achieved when the necessary, but not 
sufficient, condition 11311 



Tr(P 2(0l ) 



< 



< 



F 



rji Tr(P 1 p 2 ) 



(2.10) 



is fulfilled, where the operators Pi and P2 are the projec- 
tors onto the supports of p± and p 2 , respectively. It has 
been pointed out that there exist mixed states for which 
the fidelity bound cannot be reached for any value of the 
prior probabilities [H, [H[ . 



B. The canonical representation of the density 
operators 

When we want to explicitly determine the optimum 
detection operators, it is crucial to use convenient basis 
vectors for representing the two given states. From now 
on we focus our interest to the problem of distinguish- 
ing two states of rank d the supports of which jointly 
span a 2d-dimensional Hilbert space, because it has been 
shown that the unambiguous discrimination of two ar- 
bitrary states can be reduced to this standard problem 



[lTj| | . We start from the spectral representations for the 
two given states, 



Pi 



h\ri)(h 



1=1 



Pi 



(2.11) 

The projectors onto the supports of the states then read 

d d 

P^ElnXnl, p 2 = ElWil- ( 2 - 12 ) 



i=l 



As will become obvious later, for our purposes it is ad- 
vantageous to perform two separate unitary basis trans- 
formations in the two Hilbert spaces Til and TI2 spanned 
by the supports of p\ and p 2 , respectively, yielding two 
new sets of orthonormal basis states that are denoted by 
{|ri)} and and have the property that 



(n\sj) 



(sj\n) = Ci8, 



131 



< d < 1. 



(2.13) 
(2.14) 



Basis states of this kind have been used already previ- 
ously to study the unambiguous discrimination of two 
mixed states [ill, Ell arid to construct a very simple ex- 
ample 13]. After the basis transformations have been 
performed, the density operators take the form 



Pi 



Pi 



i,j=l 



d 

E 



\8i)( 8j \. (2.15) 



In the following we shall refer to Eqs. (|2.15p together 
with Eqs. (|2.13|) and (|2.14p as the canonical representa- 
tion of the two given density operators. 

In order to show that for any two density operators of 
rank d jointly spanning a 2d-dimensional Hilbert space 
the canonical representation always exists, and to give 
also a recipe how it can be constructed, we rely on the 
treatment given in Ref. 27]. First we observe that the 
operator P1P2P1 is Hermitean, and that its eigenstates, 
which we denote by |r^) , therefore span a complete d- 
dimensional orthonormal basis in Tli . Because P2 and Pi 
are projectors, it follows that {ri\P%\ri) — (ri|PiP2Pi|ri). 
Clearly, the norm of the state P^i) is not larger than 1, 
and moreover it is non-zero since the joint Hilbert space 
spanned by the supports of the two density operators is 
assumed to be 2d-dimensional. Hence we can establish 
the eigenvalue equation 



PP 2 Pih) =PiP>h) =C?\ri), 



(2.16) 



Now we introduce 



where < C L < 1 and (ri\rj) = 
the normalized states in Ti.2 that are given by [2j 



C 3 

and obey the equations 



CiCj 



(r i \P 1 P 2 P 1 \r j ), 



(n\sj) 



c 4 



(nlPiPsPilr, 



(2.17) 

(2.18) 
(2.19) 
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Taking into account that {ri\P\P2Pi\rj) = CfSij because 
of Eq. (|2.16[) . we immediately arrive at Eqs. (|2 . 1 3[) and 
([2~l"4|) . Thus we have shown that Eq. ([?T6| together 
with Eq. (|2.17[) provides the means for determining the 
two sets of canonical basis states { | ) } and {|si)}- Ob- 
viously, this requires the solution of a rfth-order alge- 
braic equation, resulting from the eigenvalue equation, 
Eq. (|2.16[) . We still remark that by making use of 
Piln) = Ci\ Si ) and Ptlsi) = Ciln), Eq. d^TSJ) can 
be transformed into Cf P2\ r i) — P^PiP^^i) which, with 
the help of Eq. (|2.17p . leads to the alternative eigenvalue 
equation P2Pi\si) — P\P2Pi\si) — Cf\si), as expected 
for symmetry reasons. 



C. Construction of the optimum detection 
operators 

Having obtained the canonical representation of the 
density operators to be discriminated, we are now in the 
position to make an explicit general Ansatz for the detec- 
tion operators IT and II2 that enable the unambiguous 
discrimination by satisfying Eq. (|2.2[) . For this purpose 
we define the states 



Ci 



Ci 



Si 



(2.20) 



where Si = y/l - Cf . Making use of Eqs. (|2~T3| and 
fl2~T4"l) it follows that 



(Vi\vj) = (Wi\wj) = 5 i:j 
and, most importantly, 



(Vi\sj) = {wi\r 3 ) = 0. 



(2.21) 



(2.22) 



The two joint sets of states {{|sj)}, on the one 

hand, and on the other hand, form two 

different complete orthonormal basis systems in our 2d- 
dimensional Hilbert space. Their mutual geometrical ori- 
entation is characterized by the relations 



( v i\ r j) = ( w i\ s j) = SiS^ 



(2.23) 



in addition to (vi\wj) = —CjSjj = — (ri\sj). In accor- 
dance with our earlier work [l3j we can now make the 
general Ansatz 



IT 



<y-ij\vi 



n a = £ AiKXtflil (2.24) 



which because of Eqs. (|2.15p and (|2.22p guarantees that 
P1II2 = /O2II1 = 0, as required for unambiguous discrim- 
ination. For these operators to describe a physical mea- 
surement, the coefficients and /Jy must be chosen in 
such a way that their eigenvalues, as well as the eigenval- 
ues of Ho, are nonnegative and not larger than 1. Using 
the expression IT = J2i j o>ijI\vi)(vj\I, where 



/ = ^(|r i )(r i |+K)(^|) 



(2.25) 



is the unity operator in the 2<i-dimensional Hilbert space, 
we can represent the operator Ho = / — Hi — n2 in the 
form 

d 

n = ^2 [(<% ~ a ij s i s j)\ r t)( r j\ +OijSiCj\ri){wj\ 

+ajiSjCi\wi)(rj\ + {Sij - a i3 CiCj - Pij)\wi}(wj\], 

(2.26) 



Moreover, from Eqs. (|2.24p and (|2.3p we obtain an ex- 
plicit expression for the failure probability, given by 

d 

Q = 1 - ^ SiSj(riiaijrji + mPijSji)- (2.27) 
»>i=i 

For brevity, in the rest of the paper we denote the diago- 
nal elements of the density operators and of the detection 
operators as 

Til = Ti: Sii = Si, (%ii = Q^, Pii = Pi- (2.28) 

Since £\ r i — 12i s i = 1j the conditions for the achieve- 
ment of the absolute minimum of the failure probability, 
Eqs. (|2.8p and (|2.9p . can be rewritten as 



2J (SiS.a 

d 

(SiSjPij - 6, 



> ' 13 ) o IJ 



}) = 0, (2.29) 



m 

^-F({C i ,r ij ,s ij }) = 0,(2.30) 

^72 



where the fidelity depends on the parameters that char- 
acterize the density operators in the canonical represen- 
tation, given by Eqs. (|2. 1 3(1 - (|2. 15j) . Clearly, in general 
the coefficients aij and fyj are not uniquely determined 
by these two equations alone, and a complete system of 
equations would have to be found, taking into account 
Eq. (|2.5p . However, under certain conditions Eqs. (|2.29|) 
and (|2.30p are sufficient for obtaining the optimum mea- 
surement, as we shall see in the following. In particular, 
this is the case when the canonical representation of the 
density operators is such that the expression for the fi- 
delity has a specific form, depending only on the diagonal 
elements tv and s,. 



D. An analytical solution for the optimum 
measurement 

We start by reconsidering a problem that has been re- 
cently explicitly solved with the help of a slightly different 
approach [l8| . We assume that the density operators are 
diagonal in the canonical representation, i. e. 



d 

Pi = J2n\n)(ri 

i=l 



P2 = 22a i \8 i )(a i \, (2.31) 
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where Eqs. (|2.13[) and (|2. 14[) hold for the eigenstates of 
the density operators. The fidelity is then readily calcu- 
lated from Eq. (|27]) as 



and Eqs. (|2~2"9"|) and (|2~3U|) take the form 



CAr. 



(2.32) 

0, (2.33) 
0. (2.34) 



A solution for the diagonal elements of the optimum de- 
tection operators can now be immediately read out. It is 
given by a* = a° and = /3°, where 



i 1 



(2.35) 

According to Eq. (|2.27p the failure probability Q does 
not depend on the nondiagonal elements of the detection 
operators when = rjtfy- and Sij = Sidij. We therefore 
conclude that in the optimum measurement 



1 



run 



a 



(2.36) 



since this requirement guarantees that and /3, can be 
made as large as possible while Ilo is still a positive op- 
erator, i. e. that the failure probability becomes as small 
as possible. Because of the condition on the eigenval- 
ues of the detection operators we have to require that 
< < 1. Therefore Eqs. (|2.35p only represent a 

physical solution for the optimum measurement when the 
ratio 772/771 falls within certain intervals. After replacing 
the coefficients a° and (3° outside these intervals by their 
values at the boundaries, in order to make Q as small as 
possible, we arrive at 



112. < J_ I 'Ll 

m — d * 1 ■ 



opt 

i 


= 1. 




= 


if 


opt 

i 


= «?, 




= # 


if 


opt 

i 


= 0, 


(3° pt 


= 1 


if 



Oi 



IT* < / V2 



in accordance with Ref. [if 
operators are then given by 



(2.37) 

The optimum detection 



n° pt = 



i=l 



opt 



Vi)(vi 



n 



opt = v^ pt l^)Kl, 



=1 



(2.38) 



and 



i=l 

+a° pt ^a|^)(r,| + (1 - a opt C, 2 - /?° pt )K)KIL 



(2.39) 



where in the latter expression Eqs. (|2.26p and (|2.36[) have 
been used. In order to show that these operators indeed 
describe a physical measurement, we still have to verify 
that Ilo is a positive operator. From Eq. (|2.39[) it be- 
comes obvious that Ilo can be represented by a matrix 
which consists of d decoupled two by two matrices. Tak- 



ing into account that Sfa°(3° 



[3? we find after 



minor algebra that for each of these matrices one eigen- 
value is zero and the other is given by 



A, = 



if Ci./a<./a< 1 /a, (2.40) 



or by A,; = 1 otherwise [18|. It is easy to check that the 
condition < Xi < 1 is indeed fulfilled for the eigenvalues 
Xi of the operator Ilo. 

A few direct conclusions can be drawn from the Eqs. 
(|2.37p . Obviously, when for the given prior probabilities 
of the two mixed states there does not exist a single value 
of i for which the condition in the middle line of Eq. 
(I2.37[) is fulfilled, then the optimum measurement is a von 
Neumann measurement, where the detection operators 
are projectors. In this case the failure probability of the 
optimum measurement is given by 



? op t = 1-771 Ei=i s i r * if 
? op t = 1-772 Y? l= i s ? s i if 



< Min, 



(2.41) 



In all other cases the optimum measurement is a gener- 
alized measurement, but only when the condition in the 
middle line of Eq. (|2.37j) is fulfilled for each single value 
of i, (i = 1, . . . , d), the fidelity bound of the failure prob- 
ability is obtained. Thus we have that Q pt = 2^/771 772 F 
if 



{a } < ^ < Mm {^~}. (2.42) 



Maxi 



Clearly, when Max* { C i-\f%} > M ™i { ^ } thc 
condition given by Eq. (|2.42| can never hold true and 
the overall lower bound of the failure probability cannot 
be reached. 

It is important to observe that the solution expressed 
by Eqs. (j236")l and (|2~3"?|) holds whenever the fidelity 
takes the form given by Eq. (|2.32|) . since due to Eq. 
(|2.36p the nondiagonal density matrix elements and 
do not enter the Eqs. (|2T2"9]) and (j2~3T)]) . In Sec. 
Ill we apply this solution to the optimum unambiguous 
discrimination of two particular density operators that do 
not have to be diagonal in the canonical representation. 



III. DISCRIMINATION OF STATES 
BELONGING TO A CLASS OF SIMILAR STATES 

A. The canonical representation and the fidelity 

Now we turn our attention to the unambiguous dis- 
crimination of two mixed states p\ and P2 of rank d 
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that are connected via a unitary transformation in the 
2e?-dimensional Hilbert space spanned by their joint sup- 
ports, 



f>2 



U pi u\ 



(3.1) 



where U* = U^ 1 . Since we want to determine the op- 
timum measurement by means of applying Eqs. (|2.29|) 
and (|2.30p , we first have to express the condition on the 
states within the framework of the canonical representa- 
tion. By inserting the respective density operators, given 
by Eqs. ([2TT5]) . into Eq. ([3TT]) . we obtain 



f>2 



= s ij\ s i)( s j\ = njU\r t )( r: j\U\ 



where (rAsj 



CiSij and (si|s. 



(3.2) 



The op- 



erator U transforms any state in the support of p\ into a 
state in the support of p 2 which means in particular that 
U\n) = J2k c ik\ s k), where J2k c ikC* k = (n\r 3 ) = 5 tj . In 
general, the calculation of the fidelity of these two mixed 
states is a difficult problem and cannot be performed an- 
alytically. In the following we therefore restrict ourselves 
to a special class of unitary transformations. 

We assume that the unitary transformation U can be 
decomposed into d independent unitary transformations 
Ui that act in the d mutually orthogonal two-dimensional 
subspaces spanned by the pairs of nonorthogonal states 
|r») and \si). In each of the subspaces a particular or- 
thonormal basis is given by the states |r,-) and \wi), where 



W; 



(IsJ-ain)). 



(3.3) 



Since according to Eqs. (|2 . 1 3[) and (|2.14p the inner prod- 
ucts of any two states in the combined set of the basis 
states of the two density operators are real, the class of 
transformations we consider is described by [27] 



U = U x {9 l )®U 2 {9 2 ) 



U d (9 d ), (3.4) 



where the transformations in the subspaces are rotations 
by the angle 9i, 

Ui{0i) = exp [0i (\wi){n\ - \ri){wi\)}. (3.5) 

As can be verified by expanding Ui{9i) in terms of powers 
of 9i, this is equivalent to 



_ J cos 9i\ri) + sm0i|u>i) if i = j 
r j) ifi^i- 



(3.6) 



In order to obtain the canonical representation of the 
density operators, we have to determine the eigenvalues 
and eigenstates of the operator P\P 2 P\, see Eq. (|2. 16|) . 
The projectors onto the supports of p\ and p 2 read 



where the expression for P 2 follows from the right-hand 
side of Eq. (|3.2[) . By applying Eq. (|3.6[) we easily find 
that 



P1P2P1 =^cos 2 ^|r,)(r l |, 



(3.8) 



and Eq. (|2. 16|) therefore immediately yields 

Ci = cos Oi. (3.9) 

From Eqs. ((3~H)) and Eq. (j3~5)) we then obtain 

U\n) = \ Si ) (3.10) 

which means that 



(r \U\n 



(3.11) 



After calculating the matrix element (ri\p 2 \rj) from both 
expressions in Eq. (|3.2p , using Eq. (13. 1 1|) . we finally get 

Sij = nj. (3-12) 

Hence under the condition given by Eq. (|3.4p our starting 
equation, Eq . (|3.2|) . can only be fulfilled when 



Pi=Y r v\ r i)i r jl 92 = Y r v\ s i)( s j\- ( 3 - 13 ) 



In other words, the two mixed states we consider dif- 
fer by the orientation of their respective canonical basis 
states in the 2d-dimensional Hilbert space, but the rela- 
tive weights of these states and the coherences between 
them are the same. 

After having specified the relation between the matrix 
elements of the two density operators, our next step be- 
fore applying Eqs. (|2.29p and (|2.30p is the calculation of 
the fidelity. From Eq. (|3.ip we obtain ^/p^ = U^fp~[U^ 
and Eq. (|2.7p therefore yields 

F = TrlVpTC/VpT^I 

- TrKV^riVpT^ U^pltf^pl) 1 *] 

= Tr|VpiZVpI|, (3.14) 

where we made use of the fact that WU = I. Writing 
the unity operator in our 2<i-dimensional Hilbert space 

as I = J2t=i(\ r i)( r i\ + \ w i)( w i\) an d inserting it twice, 
taking into account that = 0, we obtain 

= Tr| J2 Gi VFAn){n \ VpTI, (3.15) 



where Eq. (I3.1ip has been used. Defining the vector 
\ a i) = \/Pi\ r i), we find that F = J2i Ci Tr (l a i) ( a t\) and 
arrive at the final result 



F = Y i C i (r i \ Pl \r i ) =J2Ci* 



(3.16) 
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Interestingly, for the class of states we consider the fi- 
delity does not depend on the nondiagonal elements of 
the density operators in the canonical representation, no 
matter what is the kind of the individual unitary trans- 
formations in the two-dimensional subspaces. 



B. The optimum measurement 

We are now prepared to determine the measurement 
for the optimum unambiguous discrimination. Upon in- 
serting the expression for the fidelity, Eq. (|3.16p . into 
our basic conditions, Eqs. (|2.29p and (|2.30p . taking into 
account that Sij — , we arrive at the two equations 



SiSjf3ij r ji + 



E 



S?c 



$i Pi 




0. 



(3.17) 



that have to be fulfilled by the coefficients determining 
the optimum detection operators. Because of the special 
structure of these equations, resulting from the specific 
expression for the fidelity, we are free to make the Ansatz 



— OliSij , {3ij — Pi^i' 



(3.18) 



Obviously, the problem to be solved is then reduced to 
the problem expressed by Eqs. (|2.33[) and (|2.34[) . in the 
special case that = Sj. The previous solution, given 
by Eqs. (|2.35|) - (|2.37p , therefore can be immediately 
applied and the optimum coefficients read 



opt 

OL 

% 


= 1, 


poet 


= 


if 


opt 

OL 

% 


= <*?, 




= Pt 


if 


opt 

OL 

I 


-o, 




= i 


if 


where 












= s?( 


H 







vi — 



Ci 



< 



/2a < J_ 

rn_- Ci' 



(3.19) 



^ 1 



i 

(3.20) 

The solutions for the optimum detection operators follow 
by inserting the optimum coefficients into Eqs. (|2.38p and 

o. 

In order to obtain compact results for the minimum 
failure probability, Q op t, ensuing from the optimum mea- 
surement, it will be useful to adopt the convention that 



Ci < C 2 < . . . < C d _x < c d . 



(3.21) 



After inserting Eqs. (|3 . 1 8[) - (|3.20p into the equation for 
the failure probability Q, Eq. (|2.27p . taking into account 
that Ti = Si, we find again that the structure of the 
resulting expressions depends on the ratio of the prior 



probabilities. If the latter is such that one of the two von 
Neumann measurements is optimal, the minimum failure 
probability takes the form 



Q 



opt 



;opt 



if 

if 



2a <Ci, 
m - cr 



On the other hand, with respect to the saturation of the 
fidelity bound we find that 

Qopt = 2^j^F if C d < ^ < ±, (3.23) 

where F = Yli=i ^i r i- I n the intermediate regions of 
the ratio of the prior probabilities the optimum failure 
probability can be written as 

k d 

Qopt = i-J2( 1 - 2 VmmCi)n-m E s ? r * ( 3 - 24 ) 

i=l i=k+l 

if C k <Jf i <C k+1 (l<fc<d-l) 



(3.25) 



and 

k d 

g op t = 1-^(1-2^^-1)2 S * 



i=k+l 



if 



< 



')i — c fc V — — > 



Clearly, in dependence on the ratio of the prior probabili- 
ties of the states, there are in general 2d+l parameter re- 
gions in which the optimum measurement operators have 
a different structure and consequently the expression for 
the optimum failure probability takes a different form. 
These regions do not depend on the matrix elements of 
the density operators, but only on the canonical angles, 
that is on the constants C L . For d = 2, the calculation of 
C\ and C2 can be easily performed analytically by means 
of Eq. (|2.16p since it only amounts to the solution of a 
quadratic equation. 

It is interesting to compare the parameter interval 
in which the fidelity bound of the failure probability 
can be actually achieved, specified in Eq. (|3.23p . with 
the respective parameter interval following from a nec- 
essary, but not sufficient condition [l3j . as given in 
Eq. (|2.10p . Representing Pi and P2 as \ r i)( r i\ 
and J2i=i \ s i}( s i\i respectively, we find that Tr(P 2 pi) = 
Tr(P2/Oi) = 5Zi=i C? r i- The former interval is necessarily 
not larger than the latter, the relative difference between 
the intervals obviously being characterized by the ratio 
Si=i Cfri/ Y2i=i CiCdTi, where the explicit expression 
for the fidelity has been taken into account. 

Two special cases are worth mentioning. In the first 
one the two mixed states have equal prior probabilities 
to occur, 771 = 772 = 0.5. Since the inequality Cd < 1 < 
1/Cd certainly holds for any Cd = cos8d < 1, it becomes 
obvious from Eq. (|3.23p that in this case the fidelity 
bound of the failure probability can always be reached. 
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The second special case refers to identical canonical 
angles, Cj = cos 9 for i = l,...,d which means that 
the two density operators are connected via a rotation 
by the angle 9. We mention that for a nonorthogonal 
angle 9 this is exactly the condition that has been de- 
rived in Ref. [27| as the prerequisite for secure quantum 
communication when the two-pure-state protocol [26j is 
extended to two mixed states. In this case it follows that 
Tr(Pip2) = Tr(P2Pi) — F 2 — cos 2 9 and our general so- 
lution, represented by Eqs. (|3.22[) - (|3.25|) reduces to 

Qm= hwsF «F<Jz<i (326) 

{ Vmin + ^max-F otherwise, 

where ?7 m in(?7max) denotes the smaller (larger) of the prior 
probabilities. This result exactly corresponds to the so- 
lution for the optimum unambiguous discrimination of 
two pure states [5|. 

IV. CONCLUSIONS 

In this paper we have shown that an analytical so- 
lution for the optimum unambiguous discrimination of 
two mixed states can be obtained provided that the ex- 
pression for their fidelity is given by Eq. (|2.32|) . where 
the density operators are represented with the help of the 
canonical basis that separates the joint Hilbert space into 



d mutually orthogonal two-dimensional subspaces. The 
discrimination problem is then mathematically equiva- 
lent to distinguishing pairs of pure states. We applied 
the solution to the discrimination of two mixed states 
that belong to a special class of similar states. The den- 
sity operators of these states do not have to be diagonal 
in the canonical representation. Our results might be 
also of interest for quantum cryptography, where states 
of the kind considered in this paper play a role [27j ■ 

We still note that after finishing this work a related 
paper appeared [19( where the authors investigate lower 
bounds of the failure probability by introducing a differ- 
ent kind of state vectors spanning the joint Hilbert space. 
These states arc defined by the requirement that for any 
two density operators the expression for the fidelity is 
of a form equivalent to Eq. (|2.32[) . In contrast to the 
canonical basis states considered in the present paper, 
the states introduced in [l9| do not necessarily provide 
an orthogonal basis in the supports of the two density 
operators. 
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